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ABSTRACT 

Hydraulic  ram  concerns  the  dynamic  loading  and  cata- 
strophic failure  of  liquid  filled  fuel  tanks  impacted  by  high 
speed  projectiles.   Hydraulic  ram  is  divided  into  two  phases: 
shock  phase  and  cavity  phase.   The  shock  phase  was  studied. 
Theoretical  predictions  of  shock  radius  versus  time  were  com- 
pared with  experimental  results.   The  theory  was  found  to  be 
reasonably  accurate  for  blunt  shaped  projectiles. 

Included  in  the  theory  were  predictions  of  shock  pressure 
and  pressure  profiles  behind  the  shock.   Theoretical  pres- 
sures are  presented  graphically  but  were  not  validated  by 
experiment. 

Velocity  decay  of  the  projectiles  in  the  liquid  was  in- 
vestigated and  correlated  with  theory  to  provide  information 
for  further  studies  of  the  cavity  phase. 
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I.   INTRODUCTION 

Non-nuclear  survivability  of  aircraft  flying  in  a  hos- 
tile arena  has  been  studied  throughout  the  history  of  avia- 
tion.  Recent  experiences  in  Vietnam  have  demonstrated  the 
vulnerability  of  sophisticated  aircraft  to  both  small  arms 
ground  fire  and  the  more  complex  surface-to-air  defenses. 
Except  for  the  rare  case  of  a  direct  impact  by  a  surface- 
to-air  missile,  the  damage  has  been  incurred  by  relatively 
small  projectiles. 

Tactical  aircraft  fuel  tanks  or  cells  have  the  largest 
surface  area  and  volume  of  all  vulnerable  components.   For 
this  reason  fuel  cells  have  been  most  susceptible  to  damage. 
Projectile  impact  with  fuel  cells  have  caused  fuel  starva- 
tion, fire  or  explosion  as  well  as  extensive  structural  dam- 
age.  The  survivability  of  aircraft  when  a  ballistic  threat 
causes  destruction  of  fuel  filled  tanks  has  been  the  main 
impetus  behind  this  study.   Astronomical  costs  of  combat 
aircraft  have  made  it  cost  effective  to  acquire  more  knowl- 
edge, of  the  destructive  mechanism  to  fuel  cells.   Another 
application  of  this  information  could  be  in  space  flight 
where  meteoroid  collisions  are  possible. 

Projectiles  which  penetrate  fluid  filled  cells  cause 
damage  many  times  more  severe  than  that  incurred  by  impact 
with  an  empty  cell.   With  sufficient  projectile  energy  there 
will  be  catastrophic  failure  of  cell  structural  components. 
Hydraulic  ram  is  the  term  used  to  describe  such  destructive 
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events.   The  hydraulic  ram  phenomenon  is  usually  described 
by  two  phases:  the  shock  phase  and  cavity  phase.   The  effects 
of  these  phases  are  assumed  to  be  sequential,  occurring  in 
the  order  listed  above.   Upon  initial  projectile  impact, 
energy  is  transferred  to  the  fluid,  creating  a  strong  hemi- 
spherical shock  centered  at  the  point  of  impact.   The  shock 
pressure  loading  may  cause  failure  of  the  entrance  wall. 
Pressures  load  the  wall  for  a  very  short  time  because  of 
geometrical  attenuation.   Essentially  there  is  an  impulsive 
load  on  the  front  wall  during  the  shock  phase. 

As  the  projectile  travels  through  the  fluid,  a  cavity  is 
formed.   This  phase  has  sometimes  been  called  the  bubble 
phase  because  of  the  bubble-like  appearance  of  the  cavity. 
Others  have  called  this  phase  the  drag  phase,  using  the  la- 
bel to  describe  the  mechanism  which  creates  the  cavity  [Ref . 
1] .   Oscillation  of  the  cavity  volume  produces  pressure 
pulses  that  may  be  sufficiently  intense  to  cause  catastroph- 
ic failure  of  cell  components.   The  intensity  of  the  result- 
ing pressure  pulses  is  weaker  than  the  shock  phase  pressure 
pulse  but  the  duration  is  much  longer  [Ref.  2] .   Ultimate 
destruction  of  tank  components  other  than  the  entry  wall  has 
been  generally  attributed  to  the  cavity  phase. 

Until  recently,  experimental  data  from  studies  of  hy- 
draulic ram  have  not  been  able  to  describe  the  basic  phenom- 
enon.  Entire  aircraft  or  wing  sections  have  been  destroyed 
while  studying  the  phenomenon.   This  method  of  experimenta- 
tion has  become  too  costly;  hence  theory  explaining  the 


phenomenon  is  necessary  to  defeat  hydraulic  ram.   Much  of 
the  theory  as  it  applies  to  hydraulic  ram  is  still  in  its 
infancy.   The  complete  mechanism  is  complex,  requiring  ac- 
curate and  detailed  analytical  efforts.   Correlation  of 
theory  with  experiments  must  be  accomplished  simultaneously 
to  assess  theory  validity. 

Research  at  the  Naval  Postgraduate  School  has  been  de- 
signed to  isolate  and  observe  the  individual  aspects  of  hy- 
draulic ram.   The  ballistic  range  has  been  used  to  study 
energy  losses  through  aluminum  plates,  shock  propagation  and 
the  cavity  phase.   Specific  details  were  observed  with  the 
intent  of  applying  the  knowledge  gained  to  the  design  of  a 
survivable  fuel  cell. 

This  study  has  investigated  the  shock  phase  using  the 
theory  of  Yurchovitch  and  Lundstrom  [Ref s.  1  and  2]  .   A  com- 
puter program  was  formulated  from  the  mathematical  model. 
Experimental  data  for  various  projectile  energies  were  com- 
pared with  the  results  of  the  theory.   Yurchovitch1 s  formu- 
lation of  shock  radius  as  a  function  of  time  has  agreed  with 
the  experimental  data.   These  tests  also  have  shown  that 
damage  to  the  tank  front  wall  was  not  extensive  enough  for 
significant  fuel  leakage  at  the  energy  levels  tested. 
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II.   ANALYTICAL  STUDY 

In  order  to  predict  the  effects  of  the  shock  phase  of 
hydraulic  ram,  it  is  necessary  to  study  how  the  projectile 
impact  energy  affects  the  fluid  that  is  disturbed.   It  was 
assumed  that  the  impact  energy,  E  ,  was  known  and  was  suffi- 
cient to  create  a  strong  hemispherical  shock.   The  perturba- 
tion in  the  fluid  can  then  be  described  by  assuming  that  the 
shock  radius  is  proportional  to  time  to  the  0.8  power,  and 
that  there  is  a  power  law  density  profile  behind  the  shock. 
This  formulation  was  suggested  by  Yurchovitch  [Ref .  2] . 
Yurchovitch  found  that  for  a  given  fluid  equation  of  state 
the  energy  transferred  to  the  fluid  determines  the  resulting 
flow  field. 

A.   SHOCK  FRONT  PRESSURES 

The  one-dimensional  Rankine-Hugoniot  equation  for  a 
shock  moving  through  undisturbed  fluid  was  used  to  determine 
shock  front  conditions.   To  complete  the  solution  an  equation 
of  state  was  needed.   The  Tait  Equation  for  isothermal  com- 
pressibility of  liquids, 

-  I(3V)  =  1 ,  (2.1) 

v  8ps   n[P+B(s)  ] 

where  n  is  a  constant  for  a  specific  liquid  and  B  is  a  func- 
tion of  entropy  only,  was  integrated  along  lines  of  constant 
pressure  and  entropy.   Integrating  equation  (2.1)  yields 
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P  =  B(S)  [(v(TfO))n_1]  (2  2) 

V"(T,P) 

Assuming  the  changes  across  the  shock  are  adiabatic  and  the 
pressure  in  the  undisturbed  fluid  is  small  compared  to  the 
shock  front  pressure,  B  is  given  as 

B  =  M°    where  C*  =  &)  (2.3) 

n  °    dp  s 

With  the  assumption  that  B  and  n  are  constant  for  a  given 

fluid  the  equation  of  state  is 

Poco   P   n 
P  =  -2-£[(— )n-lJ  (2.4) 

n    pQ 

For  water  n,  a  constant,  was  found  to  be  7.0  [Ref.  3]. 

Non-dimensionalization  of  the  variables  was  found  to  be 
convenient  and  the  following  scheme  was  adopted: 

F-    r 

?"  RjTtT 

Ms  =  R8(t)/C0 

P(r,t) 


<!>(?, MJ  = 


po 


K£,MC)  =  U(r#t)  (2.5) 


Rs(t) 


f(?,HJ  =  P(f't} 


f>oHM 


g(?,Ms)  =  ??(r-t) 
s    Ri(t) 
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With  the  assumption  that  PQ<<P  and  e  <<e  the  Rankine-Hugoniot 
equations  become 

1  =  *,(H,) 

fi    =  <J>!  (2.6) 

gi  =  |f,(l-i-) 

and  the  equation  of  state  is 

f  =  l.(^n_1}  (2>7) 

nMj 

Solving  the  above  equations  for  the  parameters  at  the  shock 

the  following  is  obtained: 

■  *i  -  — r—R —  (2-8) 

i-i(^-i) 

nM| 

*i    =  -^-(^-D  <2-9> 

nMs 

fi    =   <Ci  (2.10) 

By  assuming  values  for  t\> x ,  the  equation  (2.8)  can  be  solved 
for  the  shock  Mach  number,  Ms.  Subsequently,  fj  and  <j>  j  ver- 
sus M  can  be  found.  Conditions  at  the  shock  front  are  now 
known  while  conditions  behind  the  shock  front  are  still  un- 
known. 

B.   PRESSURE  FIELD  BEHIND  THE  SHOCK 

The  density  profile  behind  the  shock  is  assumed  to  have 
the  form 

*  =  n>Ieq(Ms)  (2.ii) 
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To  find  q(M  )  the  total  mass  enclosed  by  the  shock, 
s 

p0R|2TrjWdS   , 

is  equated  to  the  original  mass  occupied  by  the  enclosed 
volume, 

2/3ttR^Po 
with  the  result 

}i^2d£  =  1/3  (2.12) 

Substituting  the  assumed  density  profile  into  (2.12) >  q(M_) 

St 

is  found  to  be 

q(Ms)  =  3(^-1)  (2.13) 

and  the  density  function  becomes 

*  =  ^^3^rl) 

It  is  now  possible  to  solve  for  the  particle  velocity,  <f>, 
by  substituting  equation  (2.12)  into  the  one  dimensional  con- 
tinuity equation  and  solving  for  the  constant  with  4*==4>  x »  at 
5=1  (the  shock  front  location) .   The  continuity  equation  is 
given  by 

(*-u3? +  *4 +  2  i  -  -ms^-s 

•  •      • 

where  3  =  RSRS/Rc  -"-s  ^he   shock  decay  coefficient.   With  the 
above  substitution  and  integration,  the  particle  velocity 
profile  is  found  to  be 

*  -  M 
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By  substituting  the  density  and  particle  velocity  pro- 
files into  the  momentum  equation 

vv  ^'H    pw        4,  3£     p  s3m 
and  integrating,  the  pressure  profile  is  found  to  be 

where  the  constant  of  integration  was  solved  with  f=f  at 
?=1. 

The  shock-  decay  coefficient,  6,  can  be  found  using  the 
assumption  that  shock  radius  is  given  by 

Rs  =  HtN  (2.14) 

Experimental  work  done  by  NASA  [Ref .  4]  has  found  that 

N  =  0.8 
for  water.   Substituting  equation  (2.14)  into  the  expression 
for  8/  the  decay  coefficient  becomes 
8  =  -0.25 
It  is  now  possible  to  solve  the  pressure  distribution 
function  given  the  shock  Mach  number,  M  . 

C.   SHOCK  MACH  NUMBER 

By  considering  the  initial  energy  being  released  at  the 

center  of  symmetry  in  the  fluid  it  is  possible  to  equate  the 

initial  energy  to  the  integral  energy  density  behind  the 

shock.   The  energy  liberated  initially,  EQ,  is  distributed 

such  that 

Rq 
EQ  =  /  (e-e0+")p4Trr2dr  (2.15) 
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Using  the  assumptions  of  a  strong  hemispherical  shock  as  well 
as  the  density  profile  and  particle  velocity  profile  as  de- 
rived previously,  the  integral  equation  (2.15)  in  non-dimen- 
sional form  becomes 


Eo  =  ^PoCoRsMs{lfr^r^?-1)*i+*?^^^q+2d5       (2.16) 

The  values  of  p0,.C0  and  n  can  be  found  for  a  given  fluid. 
By  choosing  a  Mach  number  the  dependent  variables  $l  and  \J;1 
are  known.  Equation  (2.16)  can  then  be  integrated  and 
2EQ/R|  versus  shock  Mach  number  can  be  plotted.  Given  the 
energy  released,  EQ,  it  is  then  possible  to  find  the  shock 
Mach  number  as  a  function  of  shock  radius  from  the  plot  of 
2E_/R|  versus  M_. 

There  is  however,  an  anomaly  at  a  shock  radius  of  zero. 
To  avoid  a  finite  amount  of  energy  in  zero  volume  or  an  in- 
finite energy  density,  reference  5  suggests  the  initial 
radius  be  given  as 

RQ  =  (   E°   )  x/3  (2.17) 

Poco2* 

Equation  (2.17)  was  then  used  as  the  starting  point  for  the 
solution  of  the  theoretical  pressures  behind  the  shock. 
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III.   EXPERIMENTAL  PROCEDURE 

The  basic  components  of  the  ballistic  range  used  for  this 
study  are  shown  in  figure  1.   A  down-range  view  of  the  bal- 
listic range  in  figure  2  depicts  the  chronograph  screen  ar- 
rangement and  test  tank.   The  two  spark  sources  above  the 
tank  were  used  as  light  sources  for  the  shadowgraphs.   Two 
rifles  were  used,  a  .22  caliber  and  a  Remington  .222  caliber. 
Commercial  Remington  .22  caliber  Long  Rifle  ammunition  was 
used.   The  .222  ammunition  was  hand  loaded  for  testing  so 
that  the  velocities  of  the  projectiles  could  be  adjusted  and 
would  be  known  more  precisely.   Sierra  bullets  classified  as 
spitzers  were  used.   Figure  3  displays  the  corresponding 
shapes  of  the  projectiles.   The  mass  of  the  .222  caliber  and 
.22  caliber  projectiles  are  listed  in  figure  4  along  with 
the  parameters  used  to  calculate  the  energy  of  each  configu- 
ration.  The  area  and  drag  coefficient  parameters  listed  in 
figure  4  were  experimental  values  used  to  calculate  the  ve- 
locity decay  that  was  experienced  during  projectile  passage 
through  the  test  tank. 

AVTRON,  No.  A914T333,  chronograph  screens  were  used  to 
provide  start  and  stop  pulses  to  Monsanto  101  B  counters. 
As  the  projectile  passed  through  the  screen  the  pulse  circuit 
was  energized.   The  first  screen  started  the  first  counter 
and  the  second  screen  stopped  the  first  counter  and  simultan- 
eously initiated  the  second  counter.   The  third  screen 
stopped  the  second  counter.   The  projectile  impact  velocity 

was  computed  using  these  times  and  distance  measurements. 
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Figure    3.      Projectile    Shapes    (L.toR.-.  .22 
caliber,  m=5.75xl0~3  lb,  .222  caliber,  m=6.43 
x|0"3|b,  .222  caliber,  m  =  7.85xl0"3  lb) 


Caliber 

.22  Long  Rifle 

.222  Remington 

.222  Remington 

Impact  Velocity,  Vo 
(ft/sec) 

1200 

2500 

2900 

Projectile  Mass,  m 
(lbm) 

5.75xlO"3 

6.43xl0"3 

7.85XI0"3 

Impact  Energy,  E0 
(in-lb) 

1600 

7493 

12323 

Experimental    Drag 
Coefficient,  Cq 

1.20 

0.45 

0.50 

Frontal  Area  of 
Projectile  After 
Traversing  Fluid,  A 
(in2, 

0.122 

0.098 

0.165 

Equivalent    Flat 
Plate   Area,  f  (  =  CDA) 
('in2, 

0.146 

0.044 

0.083 

Figure      4.      Projectile     Parameters 
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The  events  within  the  fuel  cell  were  timed  with  two  addi- 
tional counters.   A  fourth  chronograph  screen  provided  a 
signal  pulse  which  was  distributed  to  the  adjustable  time 
delay  unit  (ATDU)  and  the  two  remaining  counters  to  initiate 
them.   The  main  element  of  the  adjustable  time  delay  unit 
circuitry  was  a  Signetic,  N74123,  retriggerable,  monostable 
multivibrator  as  shown  in  figure  5.   Within  the  ATDU  were 
two  duplicate  delay  circuits  to  provide  a  100  volt  output 
pulse  to  each  spark  source.   The  time  delay  for  both  circuits 
in  the  unit  was  adjustable;  however,  the  duplicate  circuit 
delay  was  initiated  by  the  first  circuit's  output.   This  re- 
sulted in  the  delay  time  for  the  second  output  pulse  from 
the  ATDU  being  the  sum  of  the  delay  time  of  the  first  cir- 
cuit plus  the  delay  of  the  second  circuit.   The  delay  time 
range  for  the  first  circuit  was  500  to  1600  micro-seconds, 
while  the  second  circuit  provided  an  additional  delay  of  5 
to  15  micro-seconds.   The  100  volt  output  pulses  going  to 
the  spark  sources  were  also  used  to  stop  the  two  counters. 
The  first  output  pulse  went  to  the  spark  source  providing 
illumination  for  a  shadowgraph  in  the  horizontal  plane, 
while  the  second  pulse  triggered  the  spark  source  for  a 
shadowgraph  in  the  vertical  plane.   Each  counter  therefore 
measured  the  time  the  shadowgraphs  were  taken. 

Two  16-inch  collimating  mirrors  were  used  to  pass  light 
through  the  tank.   The  spark  sources  were  located  at  the 
focal  point  of  their  respective  mirrors.   Figure  6  demon- 
strates the  arrangement  of  the  tank  with  respect  to  the 
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spark  sources,  mirrors  and  photographic  plates.   The  side 
and  bottom  walls  of  the  tank  were  made  of  plexiglass  through 
which  very  little  light  attenuation  occurred.   One  side  wall 
and  the  bottom  wall  had  grid  lines  etched  in  the  plexiglass 
for  use  in  measuring  the  shock  radius  and  cavity  distance. 
High  speed  Polaroid  film,  Type  57,  was  used  for  the  shadow- 
graphs.  A  ground,  first  surface  mirror  was  used  to  deflect 
the  light  in  the  vertical  plane  90  degrees. 

The  tank  or  test  cell  had  front  wall  inner  dimensions  of 
17.5  inches  w.ide  and  17.0  inches  high  with  a  depth  of  17.0 
inches.   Figure  7  contains  a  view  of  the  tank.   The  main  mem- 
ber of  the  front  wall  was  one-inch  thick  aluminum  and  served 
as  a  mounting  plate  for  the  specimens  used  to  simulate  fuel 
cell  wall  material.   To  accommodate  the  test  material  the 
mounting  plate  had  a  circular  cut-out  15.0  inches  in  diameter, 
The  specimen  was  mounted  behind  the  mounting  plate  with  16 
bolts  arranged  in  a  circular  pattern.   A  front  wall  of  7075- 
T6  Aluminum,  .090  inches  thick,  with  a  prepunched  entry  hole 
is  shown  in  figure  8.   This  wall  was  used  in  the  experiments 
with  the  exception  of  the  few  tests  which  used  a  solid  front 
wall.   The  rear  wall  of  the  tank  was  identical  to  the  front 
wall. 

To  provide  support  for  the  front  and  rear  walls  aluminum 
angle  strips  were  placed  longitudinally  at  the  four  corners 
of  the  tank.   The  angle  strips  were  bolted  to  the  front  and 
rear  mounting  plates.   In  addition  to  providing  structural 
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Figure  8.   Front  Wall  Test  Specimen, 
Pre-punched  Hole 
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rigidity,  the  angle  strips  retained  the  one-inch  thick  plexi- 
glass tank  walls.   The  top  of  the  tank  was  open. 

Room  Temperature  Vulcanizing  (RTV)  rubber  sealant  was 
used  at  all  joints  to  make  the  tank  water  tight.   At  the  two 
highest  energy  levels  the  sealant  around  the  front  wall  had 
to  be  replaced  several  times.   Apparently  the  shock  pres- 
sures were  sufficient  to  deteriorate  the  sealant  bond  to  the 
plexiglass  and  aluminum.   Tape  was  used  to  prevent  water 
spillage  through  the  prepunched  hole  and  to  avoid  any  signif- 
icant energy  loss  during  penetration. 

During  initial  testing  the  tank  was  set  on  the  top  of 
the  stand.   It  was  discovered  that  although  the  tank  was  re- 
strained by  the  stand  in  longitudinal  and  lateral  motion, 
angular  tank  movement  was  created  by  the  moment  arm  and  the 
force  of  impact  at  the  center  of  the  tank.   This  invalidated 
the  shock  radius  data  and  the  tank  was  subsequently  bolted 
to  the  frame.   It  was  also  discovered  that  the  tank  and 
stand  would  move  approximately  eight  inches  during  a  test 
firing.   To  alleviate  this  problem  the  front  two  legs  of  the 
stand  were  bolted  to  the  floor.   Eventually  these  bolts 
worked  loose  and  finally  all  four  legs  were  bolted  to  the 
floor  and  150  pounds  of  sandbags  were  placed  on  the  stand. 

The  experiment  was  run  at  three  energy  levels.   One  en- 
ergy level  was  tested  before  continuing  with  the  next  higher 
energy  level.   By  varying  the  time  delay  with  the  ATDU, 
shadowgraphs  of  the  shock  radius  versus  time  were  taken.   It 
was  also  possible  to  measure  the  projectile  distance  into 
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the  tank  versus  time.   Figure  9  shows  a  pair  of  typical 
shadowgraphs.   The  impact  velocity  of  the  projectile  was 
computed  from  the  times  on  the  first  two  counters  and  the 
distance  between  chronograph  screens.   The  shadowgraph  times 
were  measured  by  the  third  and  fourth  counters. 

Acquiring  data  for  the  .22  caliber  was  found  to  be  very 
difficult  because  of  the  deviation  in  impact  velocities  from 
one  test  firing  to  another.   A  good  estimation  of  velocity 
was  required  prior  to  a  test  to  compute  projectile  time  of 
flight  from  the  fourth  chronograph  screen  to  the  tank  wall. 
That  time  plus  the  desired  time  for  photographing  the  shock 
radius  determined  the  delay  time  used  to  set  the  adjustable 
time  delay  circuit.   Since  the  .222  caliber  projectiles  were 
hand  loaded,  the  impact  velocities  were  more  predictable. 
However,  in  the  .222  caliber  tests  the  velocity  had  a  suffi- 
cient deviation  to  make  it  difficult  to  shadowgraph  a  shock 
radius  smaller  than  one  inch. 

The  data  were  reduced  by  first  calculating  the  projectile 
impact  velocity.   The  time  of  flight  of  the  projectile  from 
the  fourth  chronograph  screen  to  the  tank  was  then  calcu- 
lated.  This  time  was  subtracted  from  the  elapsed  time  of  the 
shadowgraphs  as  measured  by  the  third  and  fourth  counters. 
The  result  was  the  time  corresponding  to  the  shock  position 
from  initial  impact  at  the  front  wall.   Shock  radius,  meas- 
ured on  the  shadowgraph,  was  then  plotted  at  the  correspond- 
ing time.   Measurements  of  the  shock  radius  were  accurate  to 
±0.02  inches. 
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Figure  9.   Shadowgraphs  After  Projectile  Impact  Through 
a  Pre-punched  Front  Wall 
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Projectile  distance  into  the  tank  versus  time  data  were 
also  accumulated.   The  shadowgraph  at  24  micro-seconds  in 
figure  9  shows  that  the  projectile  was  not  distinctly  visible. 
Work  done  by  McMillen  [Ref.  6]  implies  that  the  projectile 
nose  is  at  the  tip  of  the  light  cusp  seen  at  the  leading  edge 
of  the  cavity  shadow.   Light  refraction  caused  by  density 
gradients  in  the  flow  field  forms  this  nearly  triangular  re- 
gion called  the  cusp. 
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IV.   RESULTS  AND  CONCLUSIONS 

This  research  was  conducted  to  verify  the  theoretical 
prediction  of  shock  radius  as  a  function  of  time.   In  all 
but  four  cases,  the  data  collected  utilized  projectiles  im- 
pacting the  fluid  through  a  pre-punched  hole  in  the  front 
wall  of  the  tank.   There  were  four  data  points  taken  where 
the  projectile  passed  through  a  solid  front  wall  before  com- 
ing in  contact  with  the  fluid.   This  was  done  at  the  conclu- 
sion of  the  experiment  to  see  if  there  was  a  departure  from 
the  trends  of  the  previous  data  and  to  obtain  a  qualitative 
estimate  of  the  extent  of  shock  phase  damage.   The  solid 
front  wall  data  were  taken  using  the  highest  energy  level 
tested  (12323  in-lb) . 

A.   SHOCK  RADIUS  VERSUS  TIME  -  PRE-PUNCHED  FRONT  WALL 

Figures  10,  11,  and  12  show  the  experimental  data  and 
the  theoretical  prediction  of  shock  radius  versus  time. 
Seldom  were  projectile  velocities  prior  to  impact  constant 
during  the  testing.   Changes  in  velocity  have  a  direct  in- 
fluence on  the  initial  energy.   The  data  were  plotted  and 
compared,  however,  assuming  constant  impact  energies.   It 
was  therefore  necessary  to  determine  the  shock  radius  sensi- 
tivity to  changes  in  energy.   A  plot,  using  theoretical  pre- 
dictions of  the  shock  radius  versus  impact  energy  for  a 
constant  time  (figure  13),  indicates  that  there  is  an  insig- 
nificant change  in  shock  radius  for  typical  experimental 
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variations  of  impact  energy  about  the  reference  energy.   The 
impact  energy  calculated  from  experimental  data  of  the  indi- 
vidual tests  was  within  ±200  in-lb  of  the  comparison  energy 
level.   Corresponding  changes  in  shock  radius  were  of  an 
order  of  magnitude  smaller  than  the  accuracy  of  the  radius 
measurements.   Therefore,  no  corrections  were  applied  to  the 
measured  shock  radii  for  experimental  differences  in  impact 
energy  about  the  comparison  energy  level. 

The  data  for  an  initial  impact  energy  of  1600  in-lb  had 
the  best  agreement  with  theory.   Theory  assumes  that  the 
total  projectile  energy  is  deposited  at  the  point  of  impact. 
The  shadowgraphs  indicate,  however,  that  none  of  the  projec- 
tiles came  to  rest  at  the  point  of  impact.   At  the  higher 
energy  levels  (7493  in-lb  and  12323  in-lb)  the  data  fell  be- 
low the  theoretical  results  which  implies  generation  of 
weaker  shocks  than  those  predicted  by  theory.   This  may  be 
explained  by  considering  the  projectile  as  a  piston  in  the 
early  stages  of  shock  formation.   It  would  be  expected  that 
the  rate  of  change  of  fluid  volume  displaced  by  the  piston 
determines  the  momentum  exchange  between  the  projectile  and 
fluid,  hence  the  shock  strength.   The  shape  of  the  .22  cali- 
ber projectile  is  noticeably  more  blunt  than  the  .222  cali- 
ber projectiles  (see  figure  2).   This  means  that  the  .22 
caliber  projectile  would  have  a  greater  rate  of  change  of 
fluid  volume  displaced  in  the  formation  stage  of  the  shock 
phase.   This  would  produce  conditions  which  are  more  closely 
represented  by  the  assumptions  of  the  theory.   It  is  expected 
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that  blunt  projectiles  or  warhead  fragments  produce  shocks 
that  agree  closely  with  the  theory. 

Eventually  the  shock  becomes  acoustic  because  of  geomet- 
rical attenuation.   The  theoretical  predictions  for  the  time 
at  which  the  shock  is  acoustic  were  22  micro-seconds  at  1.71 
inches  for  the  energy  level  of  1600  in-lb,  36  micro-seconds 
at  2.81  inches  for  the  energy  level  of  7493  in-lb  and  44 
micro-seconds  at  5.40  inches  for  the  energy  level  of  12323 
in-lb.   Theory  assumed  no  counterpressure  (P  <<P) .   For  this 
reason  it  is  expected  that  the  shock  will  actually  become 
acoustic  earlier  than  theory  predicts.   The  experimental  re- 
sults for  this  study  exhibited  too  much  scatter  to  accurately 
determine  when  the  shock  becomes  acoustic.   It  is  reasonable 
to  assume  the  duration  of  the  shock  phase  is  accurate  in 
order  of  magnitude.   The  cavity  begins  to  oscillate  at  a  time 
on  the  order  of  1000  micro-seconds,  which  implies  a  distinct 
time  separation  between  the  shock  and  cavity  phases.   This 
research  has  shown  that  the  shock  phase  and  cavity  phase  may 
be  studied  independently. 

B.   SHOCK  RADIUS  VERSUS  TIME  -  SOLID  FRONT  WALL 

The  four  data  points  taken  for  the  case  of  a  projectile 
passing  through  a  solid  front  wall  and  then  into  the  fluid 
are  shown  in  figure  12  as  X  symbols.   Time  limitations  pre- 
vented the  collection  of  more  data.   The  limited  results  do 
not  indicate  a  significant  departure  from  the  pre-punched 
data.   Figures  14,  15  and  16  show  typical  damage  to  the  .043 
inch  thick,  6061-T6  Aluminum  front  wall  specimen.   The  front 
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Figure  14.   Solid  Front  Wall  After  Test  (6061-T6  Aluminum, 
.043  in.  thick,  EQ  =  12323  in-lb) 


Figure  15. 


Typical  Damage  to  Solid  Front  Wall  at  Point  of 
Impact  (Front  view,  6061-T6  Aluminum,  .043 
thick,  E0  =  12323  in-lb) 


M 
C 

C 
O 

T3 
0) 
4J 

o 
s 


-p 
o 
fa 
•o 

•H 
H 
O 

en 

o 
-p 

<u 
en 

e 

Q 


VD 


0) 

Cn 

P>4 


39 


wall  experienced  plastic  deformation,  became  concave  in  the 
outward  direction  and  buckled.   There  was  some  additional 
deformation  at  a  few  of  the  restraining  bolt  locations  be- 
lieved to  have  been  caused  by  a  slightly  unsymmetrical  load 
distribution.   This  was  a  result  of  projectile  impact  not 
being  at  the  true  center  of  the  front  wall.   Projectile  pene- 
tration caused  no  visible  signs  of  cracking;  although  there 
was  permanent  deformation,  the  pressure  load  on  the  front 
wall  was  not  sufficiently  severe  to  produce  cracking. 

Soper  [Ref.  7]  calculated  the  percent  energy  change  ver- 
sus impact  velocity  for  penetration  of  various  thickness  alu- 
minum plates  with  no  fluid  foundation.   For  the  specimen  and 
impact  velocity  used  in  this  study  there  was  just  under  a  10% 
loss  of  energy  predicted  in  penetrating  the  front  wall.   Using 
the  theoretical  prediction  of  shock  radius  versus  impact  en- 
ergy in  figure  13,  there  would  be  an  insignificant  change  in 
the  shock  radius  for  such  a  10%  change  of  impact  energy. 

The  shape  of  the  projectiles  after  penetrating  a  solid 
front  wall  was  significantly  different  than  the  projectiles 
which  passed  through  a  pre-punched  hole.  This  was  the  re- 
sult of  the  copper  jacket  on  the  projectile  being  deformed 
by  both  penetration  and  the  stagnation  pressure  loading  on 
the  projectile. 

C.   PROJECTILE  VELOCITY  DECAY  IN  THE  FLUID 

The  impetus  for  investigating  the  projectile  velocity  de- 
cay in  the  fluid  is  the  requirement  to  know  the  time  rate  of 
change  of  energy  exchange  in  order  to  study  the  cavity  phase 
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[Ref.  1] .   It  is  difficult  to  predict  velocity  decay  because 
of  the  many  dependent  variables  of  the  problem.   The  analy- 
sis of  velocity  decay  suggested  by  Lundstrom  [Ref.  1]  was 
utilized.   The  equation  of  motion  used  was 

D  =  "mS  =  ^Pfv2cnA  (4'1) 

dt    2  £    D 

Assuming  a  constant  deformed  projectile  area,  constant  drag 
coefficient  and  constant  projectile  mass  it  was  possible  to 
rearrange  equation  (4.1)  and  integrate,  which  yields 

xr  -(pfCnA/2m)x 

¥_  =  e    r  U  (4.2) 

Vo 

Further  rearranging  and  integration  from  initial  conditions 
of  x=0  at  t=0  yields  the  following  expression  for  projectile 
position  as  a  function  of  time 

x  =  i-  J0n(BVot+l)  (4.3) 


where 


PfCpA 

B  "   2m 


The  density  for  the  undisturbed  fluid,  HO,  was  assumed  con- 
stant.  In  calculating  the  theoretical  results  of  equation 
(4.3)  the  value  used  for  V  was  the  projectile  velocity  im- 
mediately prior  to  impact  with  the  fluid.   It  may  be  more 
suitable  to  use  the  velocity  of  the  projectile  at  the  comple- 
tion of  shock  generation  phase  as  the  initial  velocity  in  the 
fluid,  V  .   Figure  17  is  provided  to  point  out  a  possible 
discontinuity  of  the  projectile  energy  and  velocity  profile 
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during  shock  formation.   The  initial  location  and  time  of 
this  event  would  still  be  adequately  described  by  x=0  and 
t=0.   A  value  for  the  change  in  projectile  energy  due  to 
shock  formation,  AE  ,  or  the  change  in  projectile  velocity 
due  to  shock  formation,  AV  ,  would  make  it  possible  to  solve 
for  V  .   If  the  experimental  data  had  less  scatter  than  in 
this  study  it  would  have  been  possible  to  determine  the  value 
of  AES  required  for  the  theoretical  results  to  be  coincident 
with  such  experimental  results.   It  would  then  be  possible 
to  solve  for  VQ. 

Typical  projectile  shapes  after  coming  to  rest  in  the 
fluid  are  shown  in  figure  18  for  each  energy  level  tested. 
Figure  4  shows  the  values  for  the  average  final  deformed 
area,  drag  coefficient  and  mass  used  to  calculate  the  theo- 
retical curves  in  figures  19,  20  and  21.   Holm  [Ref.  8]  has 
shown  that  the  projectiles  possess  dynamic  stability  and  do 
not  tumble,  which  contributes  to  the  deformation  experienced. 
Stagnation  pressures  on  the  order  of  lO^psi  were  found  to 
exist  on  the  projectile  surface.   These  pressures  are  large 
enough  to  cause  plastic  deformation  of  the  lead  .22  caliber 
projectiles  and  the  copper  jacketed,  .222  caliber  projectiles, 

Approximately  fifty-one  percent  (51%)  of  the  .22  caliber 
projectile  impact  velocity  was  lost  within  1.5  inches  of  the 
front  wall.   The  .222  caliber  projectiles  with  an  initial 
impact  energy  of  7493  in-lb  only  lost  17%  of  their  impact  ve- 
locity at  1.5  inches.   Projectiles  with  an  energy  12323  in-lb 
lost  25%  at  the  same  penetration  depth.   Percentage  differ- 
ences of  the  energy  lost  at  the  same  distance  from  impact 
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Figure  18.   Typical  Deformation  of  Projectiles  After 
Coming  to  Rest  in  the  Fluid  Filled  Tank  - 
Pre-punched  Hole  Results  (L.  to  R. :   .22 
caliber,  E_.  =  1600  in-lb;  .222  caliber, 


E   =  1600  in-lb; 
E   =  7493  in-lb;  .222  caliber, 
in-lb) 


=  12323 
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may  be  partially  attributed  to  the  rate  of  deformation  of 
the  projectiles.   It  is  reasonable  that  the  lead  .22  caliber 
projectile  deforms  more  rapidly  than  the  .222  caliber  jac- 
keted projectile.   It  appears  that  the  highest  energy  (12323 
in-lb)  projectile  has  reached  a  limit  at  which  the  projectile 
jacket  and  lead  core  cannot  sustain  the  enormous  stagnation 
pressure.   No  longer  was  the  deformation  similar  to  the  two 
lower  energy  projectiles  (figure  18).   The  drag  coefficient, 
deformed  projectile  area  and  percentage  of  projectile  veloc- 
ity decay  at  a  fixed  distance  for  the  two  lower  energy  levels 
appears  consistent.   At  the  highest  energy  level,  values  for 
these  variables  do  not  appear  consistent.   When  there  is  ex- 
tensive deformation,  as  in  the  case  of  the  highest  energy 
level  projectile,  the  assumptions  appear  to  be  overly  simpli- 
fied. 

The  concept  of  an  equivalent  flat  plate  area,  f=C  A,  may 
be  more  advantageous  than  considering  the  drag  coefficient 
and  deformed  projectile  area  separately.  Also  included  in 
figure  4  is  the  value  for  the  equivalent  flat  plate  area,  f, 
which  best  matches  theory  to  experiment.  This  parameter  re- 
duces the  number  of  dependent  variables  by  one  and  describes 
the  projectile  drag  characteristics  with  one  number. 

In  reference  7,  minimum  impact  velocities  to  penetrate 
various  thickness  plates  are  given.   The  experimental  veloc- 
ity decay  data  were  extrapolated  by  theory  to  predict  pro- 
jectile velocities  at  the  rear  wall  of  the  tank.   These  were 
compared  with  the  minimum  impact  velocity  to  penetrate  the 
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rear  wall.   After  traveling  the  depth  of  the  tank,  17  inches, 
the  .22  caliber  projectile  (E  =  1600  in-lb)  was  predicted 
to  have  a  velocity  of  0.394  ft/sec.   Experimentally,  the  .22 
caliber  projectiles  were  found  approximately  six  inches  from 
the  front  wall  where  equation  (4.2)  predicts  a  velocity  of 
71  ft/sec.   The  .222  caliber  projectiles  with  an  energy  of 
7493  in-lb  were  predicted  to  have  a  velocity  at  the  rear  wall 
of  290  ft/sec,  which  was  32%  of  the  minimum  impact  velocity 
required  for  penetration.   The  experiments  with  the  .222  cal- 
iber projectiles  at  an  energy  of  7493  in-lb  were  the  only 
test  cases  which  made  a  noticeable  dent  but  no  penetration 
of  the  rear  wall.   The  highest  energy  level  projectiles  were 
predicted  to  have  a  velocity  of  105  ft/sec  and  there  were  no 
signs  of  projectile  impact  with  the  rear  wall.   These  results 
seem  to  indicate  that  the  exponential  velocity  decay  is  rea- 
sonable where  projectile  deformation  is  not  extensive. 

D.   PRESSURE  VERSUS  RADIUS  AND  TIME 

Figures  22,  23,  24  and  25  show  theoretical  predictions 
of  the  shock  phase  pressure  pulse  and  the  pressure  profiles 
behind  the  shock  at  various  times.   Figure  22  is  provided  to 
demonstrate  the  order  of  magnitude  of  the  shock  pressure  from 
initial  impact  until  the  shock  becomes  acoustic  for  an  energy 
level  of  1600  in-lb.   The  remainder  of  the  pressure  versus 
radius  curves  were  limited  to  a  range  of  5  micro-seconds  to 
30  micro-seconds.   This  information  was  used  to  estimate 
pressure  conditions  which  will  be  encountered  when  further 
experimental  studies  are  undertaken  to  verify  the  theory. 
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The  pressure  profiles  behind  the  shock  were  not  plotted  at 
negative  pressures  although  negative  pressures  behind  the 
shock  front  were  predicted  in  the  theoretical  results.   Since 
water  is  nearly  incompressible,  cavitation  would  occur  as  the 
pressures  became  negative.   Cavitation  has  been  observed  on 
film  records  of  similar  hydraulic  ram  studies  done  at  the 
Naval  Weapons  Test  Center,  China  Lake  and  provided  for  the 
author  by  Lundstrom. 

Immersion  of  pressure  transducers  in  the  fluid  is  being 
considered  as  a  means  of  measuring  pressures.   Figures  26, 
27  and  28  were  plotted  to  estimate  the  magnitude  of  the  pres- 
sures experienced  at  various  distances  from  the  front  wall 
and  the  duration  of  the  pressure  pulse.   The  plots  indicate 
the  need  for  a  transducer  with  a  very  fast  response  time. 

It  is  expected  that  actual  pressures  will  be  somewhat 
lower  than  predicted,  with  the  possible  exception  of  blunt 
nose  projectile  data.   There  are  two  reasons  for  predicting 
lower  pressure  profiles.   As  explained  earlier,  the  shock 
radius  versus  time  data  implied  a  weaker  shock  than  the  theo- 
retical predictions.   This  was  the  result  of  the  assumption 
in  the  theory  that  total  impact  energy  was  deposited  in  the 
fluid  during  shock  formation.   A  weaker  shock  would  result 
in  correspondingly  weaker  shock  pressures  throughout  the 
shock  phase.   The  second  reason  pertains  to  the  pressure  pro- 
file as  the  shock  Mach  number  approaches  one.   At  lower  shock 
Mach  numbers,  the  assumptions  of  undisturbed  fluid  pressure 
in  front  of  the  shock  being  much  less  than  at  the  shock 
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CPQ<<P}  will  become  invalid.  This  assumption  should  cause 
the  theoretical  predictions  to  be  higher  than  experimental 
results. 
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V.   RECOMMENDATIONS 

1.  The  scatter  of  the  experimental  data  measuring  shock 
radius  versus  time  could  be  reduced  by  placing  the  fourth 
chronograph  screen  closer  to  the  tank.   New  counters  with 
high  frequency,  temperature  compensated  oscillators  to 
provide  a  better  time  base  are  also  recommended. 

2.  The  ballistic  range  should  be  modified  to  provide  the 
capability  of  firing  higher  energy  projectiles.   Shock 
radius  versus  time  and  projectile  distance  traveled  in 
the  fluid  versus  time  should  be  studied  at  these  higher 
energy  levels. 

3.  Fluid  pressure  at  and  behind  the  shock  should  be  com- 
pared with  theory  at  various  energy  levels. 

4.  This  study  used  lead  or  copper  jacketed,  lead  core 
projectiles  and  no  tumbling  was  observed.   Further  stud- 
ies using  steel  jacketed  projectiles  should  be  completed 
to  test  for  deviations  from  the  data  of  this  study.   This 
information  would  provide  the  fuel  cell  designer  with 
data  typical  of  ballistic  threats  encountered  in  combat. 

5.  Strain  gauges  should  be  placed  on  the  front  wall  spec- 
imens to  determine  the  stresses  incurred  during  the  shock 
phase.   Strain  data  could  then  be  compared  with  SATANS 
computer  code. 

6.  Higher  energy  projectiles  which  produce  cracking  or 
catastrophic  failure  of  the  front  wall  specimen  should 
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be  tested.   These  data  would  provide  the  designer  with 
the  critical  energy  level  at  which  shock  phase  damage 
may  result  in  an  aircraft  loss. 
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COMPUTER  PROGRAM 


C  ********************************************************* 

c  *  * 

c  ♦  * 

C  *  THIS  PRCGRAM  CALCULATES  THE  FOLLOWING                 * 

C  *  PART  J                                            * 

C  *  1.  SHOCK  RADIUS  VERSES  TIME                    * 

C  *  2.  SHOCK  MACH  NUMBER  VERSES  TIME                * 

C  *  PART  II                                           * 

C  *  1.  PRESSURE  VERSE  TIN<E  AND  RADIUS               * 

C  *  FROM  IMPACT  BEhINC  THE  SHOCK                 * 

C  *  FCR  A  GIVEN  FLUID  AND  IMPACT  ENERGY, EC.   THIS  PRCGRAM* 

C  *  ASSUMES:  SHOCK  RADIUS  IS  PROPORTIONAL  TO  TIME  TO  THE  * 

C  *  C.S  POWER  LNTIL  ACOUSTIC,  A  POWER  LAW  DENSITY  PROFILE* 

C  *  BEHIND  THE  SHOCK,  A    STRONG  SHOCK,  ADIAaATIC  ACROSS    * 

C  *  THE  SHOCK,  THE  UNDISTURBED  FLUID  FBESSURE  IS  ESSEN-   * 

C  *  TIALLY  ZERO,  SPECIFIC  INTERNAL  ENERGY  IN  LNCISTURBEC  * 

C  *  FLUID   IS  ESSENTIALLY  ZERO.                           * 

C  *  * 

C  =*  * 

C 
C 

C      SYMBOL  DESIGNATION 

C  EO=  IMPACT  ENERGY  OF  PROJECTILE  (IN-LEJ 

C  C0  =  SPEED  OF  SOUND  FCR  FLUID  (FPS) 

C  RHOO=  DENSITY  IN  UNDISTURBEC  FLUID  (LE/FT3) 

C  DC=  INCREMENT  OF  NON-DIMENSIONAL  RADIUS,  DELTA 

C  ZHETA 

C  DELT  =  INCREMENT  OF  TIME  (MICRC-SEC) 

C  ET=  ENERGY/VOL.  =  2*E0/(R**3) 

C  Fl=  NON-DIMENSICNAL  PRESSURE  AT  THE  SFCCK  VS. 

C  SHOCK  MACH  NO. 

C  DODM  =  DERIVATIVE  OF  NON-DIMENSIONAL  PARTICLE 

C  VELOCITY  W.R.T.  SHOCK  fACF  NO.  VS  SFCCK 

C  MACH  NO. 

C  VS  =  SHOCK  MACH  NC. 

C  X=  PERCENT  EG  LOST  AT  ENTRY  WALL 

C 

c 

DIMENSION   ET(20) ,PSI1(20) ,F1(20) ,VS(20) ,  CCDM20) 
1,  RXMS(20) ,CRS(20J, ARMS (20) , DT(20),T(20),RE(20), 
lRVC20J-f  RDOD(20)  fXMS(20)  fRS(20)  f  PPSK20J  f  RFC20) 

DIMENSION   TT(300) ,XM(300) ,  R ( 30C ) , F2( 300 ) , FS 12 (3 00 1 , 
1DCDK2(3G0)  ,Q(300)  ,C(100)  ,F(2   ,20  )  ,P(2   ,20) 
DIMENSION  CONRSOOC)  ,  SHOO  V(  300)  ,TCF0R(3  00) 
DIMENSION  F0R(20) 
DC  500  J=1,2G 
F0R(J)=C.O 
500  CONTINUE 
N  =  18 

READ  (5,100)  EO,CO,RHCO,X,DC,DELT 
READ  (5,101) (ET( I  )  ,1  =  1, N  ) 
READ  (5,101)  (PSIim  •  1  =  1  «  N) 
READ  (5, 101) ( Fl(  I  i  ,  1  =  1, N  ) 
READ  (5,101) ( VS( I)  ,  1=1, N  ) 
READ  (5,101)(D0DM(  I)  ,  I = 1 , N ) 
WRITE  (6,200  )  EO,CO,RhGO, X, DC, CELT 
DC  15  1=1, N 
15  WRITE  (6,201)   ET  (  I  ) ,  PS  1 1  (  I  ) ,  F  1  (  I  )  ,  VS  (  I  )  ,  DCCf*  (  I ) 
C 

C      PART  I 
C 
C 
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C      PERCENTAGE  OF  INITIAL  ENERGY  DEPOSITED  ANC  INITIAL 

C      ShGCK  RADIUS 

C 

EO=X*EO 

RG=(E0*32-.2*144./(RH00*C0*C0*2.*3.1416)  )**0. 33333 
C 

C      INITIAL  VALUE  OF  ENERGY/VOL. 
C 

ET0=2.*EC/(RC**3.0) 
C 

C      INVERT  INPUT  ARRAYS  FOR  USE  WITH  INTERPOLATION 
C      ROUTINE,  PIF2 
C 

DC  12  1=1, N 

J=N-I+1 

PE(I)=ET(J) 

RPSKI  )=PSI1(J) 

RF(I)=F1(J) 

RVtI)=VS(J) 
12  RDGD(I)=DGDM(J) 
C 

C      INTERPOLATE  FOR  INITIAL  VALUES  AT  SHOCK 
C 

VS0=PIF2(ET0,RE,N,RV) 

F10=PIF2(VSO,RV,N,RF) 

PSI10=PIF2(VS0,RV,N,RPSI) 

DCDK0=PIF2 (VSO,RV,N,RDCC) 
C 

C  COMPUTE     INITIAL    TIME    FOR    DISTANCE    RO 

C 

TC=R0/(CC*VS0*12.0) 
C 

C      ASSIGN  SUBSEQUENT  VALUES  OF  SHOCK  MACH  NO.  /ST  WHICH 
C      SHCCK  RADIUS  WILL  BE  COMPUTED 
C 

DC  10   1=1, N 

IF  (VSO-VS(T))   1,2,2 
1  GC  TO  10 
I    M=N-I+1 

XMS(1)=VS0 

GC  TO  11 
1C  CONTINUE 
C 

C      COMPUTE  SUBSEQUENT  VALUES  OF  SHOCK  RADIUS 
C 

11  RS(1)=R0 

DO  20   J=2,M 

XMS(J)=VS(I) 

RS(J)=(2.*E0  /(ET(I)J  )**0.3333 

1=1  +  1 
20  CONTINUE 

K  =  M-1 
C 

C      COMPUTE  DELTA  TIME  CORRESPONDING  TO  SHOCK  MACH  NO. 
C      ANC  SHOCK  RADIUS 
C 

DC  30   1=1, K 

DRS(I)=  RS(I+1)-RS(I) 

ARMS (I )=2„/(XMS( I)+XMS(I+1H 
30    DT(I)=DRSm*ARMS(  I  )/(C0*12.0) 

T(1)=T0 

DC   40   J=2,M 
40,  T(J)=  T(J-1J+DT(  J-l) 

WRITE  (6,102) 

WRITE  (6,103)   EO 

WRITE  (6,110) 
C 

C  ,     PRINT  SHOCK  MACH  NC.  AND  SHCCK  RADIUS  VS  TINE 
C 


DC  50   1=1, M 
5C  k\RITE  (6,111)   T(  I  ),XMS(I)  ,  RS  (  I  ) 
100  FORMAT   (6E10.0) 
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1C1  FORMAT   (8E10.0) 

102  FORMAT  (1H,//,39X,  '  ENERGY ' , 15X,  'FLUID'  ) 

103  FORMAT  (29X,1E20.7) 

110  FORMAT  (1H,//,39X  • TI  ME  '  , 17X , • MACF •  » 17X  ,  • FACIUS •  ) 

111  FORMAT   (29X.3E20.7J 
2CC  FORMAT  (//,  &E20.7) 
2C1  FORMAT  (15X,5E20.7) 

C 

C 

C      PART  II 

C 

C 

C      INITIALIZE  VARIABLES  AND  ASSIGN  VALUES  USEC  WITH 

C      NON-DIMENSIONAL  PRESSURE  DISTRIBUTION  FUNCTION 

C 

BeB=T(M)*l.E06+1.05 

L  =  EBB 

TT(1)=0.0 

XM(1)=XMS(1) 

R(1)=R0 

F2(1)=F10 

PSI2(1)=PSI10 

DCDM2(1)=DCDM0 

Q(1)=3.*(PSI10-1.0) 

E=-0.25 

I=1.0/DC  +1.05 

C(1)=0.C 

WRITE  (6,301) 

WRITE  (6,302)  TT(1),XM(1) ,R(1) 

WRITE  (6,303) 

DO  60    J=2,N 
6C  C(J)=C(  J-U+OC 
C 

C      CCMFUTE  NON-DIMENSIONAL  AND  CIMENSICNAL  PRESSURE  FOR 
C      INITIAL  ShCCK  RADIUS. 
C 

DC  70  J=1,I 

F(l,J)=F2(l)+PSI2(l)*(C(J)**(Q(l)+2.0)-l.C)*(F2(l)/ 
1PSI2(1)-6*(F2(1)+XM(1)*CGDM2(1)  )  )  / (Q( 1 ) +2  .0  ) 

PU,J)=F(1,J)*RHOO*CO*CO*XM(1)*XM1)/(32.2*144.) 

PUNCH  91,R(1),C(J) ,P(1,J) 
70  CONTINUE 

K  =  l 

JJ  =  1 

SLM=0.0 

CALL  FORCECFQRtCtP.fRf  ItJJfKJ 

DC  501  J=1,I 

SLM=SUM+FGR( J) 

WRITE (6, 304)  C( J ) , P ( 1, J ) , F( 1, J ) , FCR( J) 
501  CCNTINUE 

TCFOR( 1)=SUM 
C 

C      ITERATION  TO  SOLVE  PRESSURE,  SHCCK  RADIUS,  SHCCK  MACF 
C      NO.  VS  TIME 
C 

DC  80   K=2,L 

TT(K)=TT(K-1)+DELT 

IF  (TT(K)-TO)    21,21,22 

21  XM(K)=XM(JL) 
R(K)=R(1) 
F2(K)=F2(1) 
PSI2(K)=PSI2(1) 
CCCM2(K)=D0DM2(1) 
G(K)=  Q(l) 

CO  55   J=1,I 
55  P(2,J)=P(1,J) 
PUNCH  92,R(K) 
GC  TO  80 

22  X=TT(K) 

c 

C      INTERPOLATION  OF  INPUT  DATA  AND  PART  I  DATA  FOR 
C      VARIABLES  IN  PRESSURE  FUNCTION 
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XM(K)=PIF2(X,T,M,XNS) 

R<K)=PIF2(X,T,M,RS) 

Y=XM(K) 

F2(K)=PIF2(Y,RV,N,PF) 
FSI2(K)=PIF2(Y,RV,N,RPSI) 
DCDM2(K)=PIF2(Y.RV,N,P0CD) 
Q(K)=3.*(PSI2(K)-1.C) 
WRITE   (6,301) 

WRITE   (6,302)  TT(K),XM(KJ ,R(K) 
WRITE   (6,303) 
DC  90   J  =  l  ,  I 

F(2,J)=F2(K)+PSI2(K)*(C(J)**(Q(K)+2.0)-1.0)=MF2(K)/ 
1PSI2(K)-B#(F2(K)+XM(K)*D0DM2(K) ))/(Q(K 1+2.0) 
P(2,J)=F(2,J)*RH0C*C0*C0*XM(K)*XN(K)/(32.2*144.) 
PUNCH  91,R(K),CUJ  ,P(2,J) 

90  CONTINUE 
JJ  =  2 
SUM=0.0 

CALL  FCRCE(FORfC,P,R,I,JJ,K) 
DC  502  J=1,I 
SLM=SUM+FGR(J) 

WRITE (6,304)  C(J),P(2,J),F(2,J),FCR(J) 
502  CONTINUE 

TCFCR(K)=SLM 
8C  CONTINUE 

WFITE(6,305) 

CALL  CONVER(R,XM,TCFOR,TT,L) 

91  FCRMAT  (3E15.5) 

92  FCRMAT  (1E15.5) 

201  FCRMAT  ( 1H ,//  ,36X, ' TI NE' t17Xf ■ MACF*  ,17X, • RADIUS'  ) 
302  FCRMAT  (29X.3E20.7) 

3  02  FCRMAT ( 1H , // ,39X, • P/RS • , 17X  ,  ■ P '  ,  17X  ,  « F' , 19 X, 'FORCE'  ) 
304  FCRMAT(29X,4E20.7) 

3  05  FORMAT (lh,//,36X, • T I  ME • , 17X , 'VELC « , 17X, • RACLI S • , 17X, «T 
1GF0RCE'  ) 
ENC 

FUNCTION  PIF2  ( X , XLIST, N, FL 1ST ) 
C 
C      SECOND  ORDER  INTERPOLATION 


C 


DIMENSION  XLIST  (100),  FLIST  (100) 

ELIF  (P,Q,R,S,T)  =  ( (Q-P)*(S-T)/(R-Q)+S) 

IF  (X-XLIST(N))  2,1,1 

1  I  =  N-l 
GC  TO  5 

2  IF(X-XLISTd)  )  4,4,6 

4  1  =  1 

5  K  =  1 

GC  TO  30 

6  K  =  2 

7  DO  8  I  =  1,N 

IF  (X-XLIST(I))  9,9,8 

8  CONTINUE 
I=N 

9  1  =  1-1 

30  BLIF1  =  BLIF(X,XLIST(I) ,XLIST( 1+1 ) , FL 1ST ( I ) ,FL 1ST ( 1+ 1) 
1) 

10  IF  (K-l)  11,11,12 

11  PIF2  =  BLIF1 
RETURN 

12  IF((I+2)-N)  13,13,16 

13  IF  ((I-l)-l)  15,14,14 

14  IF(ABS( XL  I  ST ( I- 1 ) -X J-ABS ( XL  I  ST ( I +  2 )-X) ) 16,15,15 

15  L  =  1+2 
GO  TO  17 

16  L  =  1-1 

17  BLIF2  =  BLIF  (X, XL  1ST ( I ) , XL  1ST ( L ) , FL 1ST (  I ) , FL I  ST ( L ) ) 
P1F2  =  dLIF  (XiXLIST( 1+1) tXLlST(L) tBLIFl , BLIF2) 

18  RETURN 
ENC 
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SUBROUTINE  FORCE ( FOR,C  ,P  ,R , I,JJ,KJ 

Q         ********************************************************* 

c  *  * 

c  *  * 

C  *  TMS  SUBROUTINE  COMPUTES  ThE  FORCE  ON  THE  ENTRY  WALL  * 

C  *  OLE  TO  SHCCK  PHASE  FOR  A  FRCNT  WALL  WITH  A            * 

C  *  PRE-PUNCHEC  HOLE, I.E.  PROJECTILE  CCES  NOT  IMPINGE  CN  * 

C  *  FRCNT  WALL  UPON  ENTRY.   COMPUTES  FCRCE  TC  BE  ZERO  FCR* 

C  *  PRESSURE  FIELD  OVER  THE  PRE-PUNCHEC  HOLE  CR  WHERE  TEE* 

C  *  PRESSURE  GOES  NEGATIVE.  THE  PROGRAM  ACCGf PL  I  SHES  THIS* 

C  *  dY  A  STAIGHT  LINE  APPROXIMATION  BETWEEN  PRESSURE      * 

C  *  DATA  POINTS.                                            * 

C  *  RHRS=  NCN-CIMENSICi\AL  RADIUS  OF  FRE-PUNCFEC  HCLE      * 

C  *  S=  SLOPE  OF  STRAIGHT  LINE  BETWEEN  PRESSURE  CATA       * 

C  *  PCINTS                                            * 

C  *  RX  =  THE  INTERCEPT  CF  R/RS , NON-C I V  ENS  ION AL  PACIUS,     * 

C  *  WHERE  PRESSURE  IS  ZERO                            * 

C  *  FCR(N)=  FORCE(LB)  FCR  AN  ANNULUS  WHOSE  OLTSIDE  RADIUS* 

C  *  IS  C(N)  AND  INSIDE  RADIUS  IS  C(N-l)           * 

C  *  * 

C   *  * 

c   ********************************************************* 

C 

c 

DIMENSION  FOR (20) , C (100) fP( 2,20) ,R(300) 
BhRS=.375/R(K) 
IF(RHRS-1. )  100,25,25 
IOC  IF(JJ-l)  1,1,2 

1  CC  10  N=2,  I 
S=(P(1,N)-P( 1,N-1) )*4. 
IF(P(1,N-D)  3,4,4 

3  IF(PdrM)  5,5,6 

5  FCR(N>=C.O 
GC  TC  10 

6  RX=(S*C(N-1)-P(1,N-1))/S 
IF(RHRS.GT.C(N) )  GC  TO  43 
IF(RX.GT.RHRS)  GO  TC  60 
PH=S*(RHRS-C(N-1) )+P(l,N-l) 
FCR(N)={.5*(P(1,N)-PF)+PH)*3.1416*(C(N)*C(M-RHRS 

1*RHRS)*R(K)*R(K) 
GC  TO  1C 
60  FCR(N)=P(1,N)*.5*3.1416*(C(N)*C(N)-RX*RX)*R(K)*R(K) 
GC  TO  1G 

4  IF(RHRS-C(N-1) )  40,40,41 

4C  FCR (N)  =  (.5*(P (It N)-P( 1,N-1)  )+P(l,N-l ))*3. 1416* (C(N) 
1*C(N)-C(N-1)*C(N-1) )*R(K)*R(K) 
GC  TO  10 

41  IF(RHRS-CCN) J  42,43,43 

42  PF=S*(RhRS-C(N-l) )+P( 1,N-1) 
FCR(N)=(.5*(P(1,N)-FH)+PH)*3.1416*(C(N)*C(N)-RHRS 

1*RHRS)*R(K)*R(K) 
GC  TO  10 

43  FCR(N)=C.O 
1C  CCNTINUE 

RETURN 

2  DC  20  N=2,  I 
S=(P(2,N)-P(2,N-1)  )*4. 
1F(P(2,N-D)  13,14,14 

13  IF(P(2,N))  15,15,16 

15  FCR(N)=G.0 
GC  TO  20 

16  RX=(S*C(N-1)-P(2,N-1) )/S 
IFtRHRS.GT.C(N))  GC  TC  143 
IF(RX.GT.RHRS)  GO  TO  160 
Ph=S*(RhRS-C(N-l) )+P(2,N-l) 
FCR(N)=(.5*(P{2,N)-FH+FH)*3.1416*(C(N)*C(N)-PHRS 

1*RFRS)*R(K)*RIK) 
GC  TO  20 
16C  FCR(N)=P(2,N)*.5*3.1416*(C(N)*C(M-RX*RX)*R (K)*R(K) 
GC  TO  20 

14  IF(PHRS-C(N-1) )  140,140,141 


65 


140  FCR(N)=(.5*<P<2,N)-P(2,N-1))+P(2,N-1))*3.1416*(C<N) 
1#C(N)-C(N-1)*C(N-1) )*R(K)*R(K) 

GC  TO  2C 

141  IF(RHRS-C(N) )  142,143,143 

142  PH=S*(RHRS-C(N-1) Wp(2,N-1) 

FCR(N)=( .5*(P(2, N)-PH)+PH)*3.141d*(C(N)*C(N)-RHRS 
l*RhRS)*P(K)*R{K) 
GO  TO  2G 

143  FCR(N)=0.0 
2C  CCNTINUE 
25  RETURN 

END 

SIEROUTINE  CGNVER ( R ,XM ,TOFGR,TT , L ) 
C 

C      TFIS  SUBROUTINE  CONVERTS  VELOCITY  AND  RADIUS 
C      TO  KM/SEC  AND  CM  RESPECTIVELY 


C 


DIMENSION  Rt3GO},XM300),TOFOR(300),CONRS(3C0),ShOCV( 
13CC),TT(30C) 

DC  510  11=1, L 

CCNRS(IIJ=R(II)*2.54 

ShOCVt II)=XM( II)*4.9/3.281 

WRITE (6,520  J  TT (  1 1  ) , SHGCV ( I  I ) , CONRS ( 1 1 ) , TCFCR ( 1 1 ) 
52C  F0PMAT(29X,4E20.7) 
510  CONTINUE 

RETURN 

END 
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